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Abstract 

We consider electroweak singlet dark matter with a mass comparable to the Higgs mass. The 
singlet is assumed to couple to standard matter through a perturbative coupling to the Higgs 
particle. The annihilation of a singlet in the mass range mj ~ mj, is dominated by proximity to 
the W, Z and Higgs peaks in the annihilation cross section. We find that the continuous photon 
spectrum from annihilation of perturbatively coupled singlets in the galactic halo can reach a 
level of several per mil of the EGRET diffuse 7 ray flux. 
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1 Introduction 

Revealing the distribution and nature of 
dark matter is one of the most interesting 
current challenges of particle physics and 
astrophysics. Numerous candidates have 
been proposed and studied, including ax- 
ions [1], neutralinos [2], Kaluza-Klein pho- 
tons [3], Kaluza-Klein or string dilatons 
[4], and superheavy dark matter either di- 
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rectly from inflationary expansion [5,6] or 
from preheating after inflation [7]. Unsta- 
ble particles like axions and dilatons have 
to be very light, in the sub GeV range, 
to survive long enough to serve as dark 
matter. However, neutralinos are usually 
assumed to have masses beyond 100 GeV, 
although lower mass limits strongly de- 
pend on supersymmetric models [8,9,10]. 

Another very interesting model for dark 
matter which can have a large mass is a sta- 
ble electroweak singlet S which couples to 
standard model matter exclusively through 
a coupling to the Higgs boson H, 
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We focus in particular on the Z 2 symmetric 
model proposed in [11,12,13,14,15], which 
also allows for a bare mass term for the 
singlet, but does not include Higgs-singlet 
mixing terms. The Lagrangian in the scalar 
sector is 
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The derivatives are the appropriate 
SU{2) x U{\) covariant derivatives acting 
on the Higgs field. We assume perturbative 
singlet-Higgs coupling for our calculations 
of singlet annihilation cross sections in the 
non-relativistic limit. We report results in 
particular for rf = 0.1 and for rf = 0.01. 
The singlet sector may also include a Z 2 
symmetric singlet self-interaction ~ X$S 4 
if the positive coupling As is weak enough 
such that its loop contributions can be ne- 
glected in the present perturbative calcu- 
lation of singlet annihilation cross sections. 
The assumption of perturbative couplings 
in the non-relativistic limit is compatible 
with the fact that the f3 functions are pos- 
itive in leading order in the couplings [14]. 
We assume that the singlet vacuum expec- 
tation value vanishes, v s = (S) = 0. Other- 
wise the singlet-Higgs coupling would yield 
a singlet-Higgs mixing term ~ r]v s Vhsh. 

Later on we will also allow for a set of 
N singlet states with a global 0(N) sym- 
metry and vanishing vacuum expectation 
values. This ensures mass degeneracy and 
universality of the singlet-Higgs coupling 
strength rj. 

The Higgs- channel between dark mat- 
ter and the standard model was denoted 
as a Higgs portal in [15]. The model pro- 
vides a minimal renormalizable dark mat- 



ter modeQ [11,12,13,14,15,16,17,18]. It 
has also been discussed as a model for 
quintessence in [22]. McDonald introduced 
the variant with a complex singlet [12], 
and effects of gauging the ensuing hidden 
U(l) symmetry are discussed in [23]. The 
presence of the singlet coupling obviously 
modifies the Higgs effective potential and 
impacts electroweak symmetry breaking 
[24], eventually triggering a strong first 
order phase transition [25,26]. 

A light electroweak singlet could reveal 
itself as missing energy in collider based 
Higgs search experiments [27,28,29,23] 
or in B decays [17]. March- Russell et al. 
pointed out that electroweak singlets and 
their fermionic partners in supersymmet- 
ric theories can also be very heavy, with 
masses up to 30 TeV [30]. 

In the present paper we will be con- 
cerned with the prospects of observation 
of intermediate mass electroweak singlets 
through their annihilation products in cos- 
mic rays. We will take the proposals of 
an electroweak singlet coupling through 
the Higgs portal minimal renormal- 
izable dark matter model seriously, and 
discuss possible annihilation signals under 
the assumption ms ~ m^. The relevant 
mechanism for a sizable signal for an in- 
termediate mass singlet would be annihila- 
tion through an intermediate Higgs boson 
into W and Z bosons. Because of the im- 
portance of the opening of the W channel 
at 80 GeV, we will denote the mass range 
80 GeV < ms < 1 TeV as the intermediate 
mass range for electroweak singlets. Prox- 
imity of m s to the W, Z, and Higgs peaks 



1 Other recently proposed classes of mini- 
mal dark matter models introduce heavy elec- 
troweak multiplets with a lightest neutral 
component [19,20], or another scalar messen- 
ger between dark matter and standard matter 
[21]. 
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in the annihilation cross section increases 
the product va substantially, thus poten- 
tially yielding a strongly enhanced flux 
of annihilation products from the galactic 
halo even without invoking boost factors 
from strong local dark matter overdensi- 
ties. However, we will see that the effect 
of the enhanced annihilation cross section 
is partially compensated for in standard 
Lee-Weinberg theory for the creation of 
thermal relics, because the requirement 
Qs = Qdm together with mj ~ will re- 
quire iV-plets of electroweak singlet states, 
and the net effect is a scaling of the flux 
j oc vcr/N. 

The strength of a direct galactic dark 
matter annihilation signal depends on the 
dark matter distribution in our galactic 
halo. The cosmic ray flux from a Navarro- 
Frenk- White (NFW) dark matter halo and 
a cored isothermal halo will be reviewed in 
section 2. Model independent limits on the 
flux will be discussed in 3. 

Annihilation cross sections which follow 
from the singlet-Higgs coupling (1) are re- 
ported in section 4. The application of Lee- 
Weinberg theory for intermediate mass 
electroweak singlets is discussed in section 
5. Section 6 summarizes our conclusions. 

2 The flux from the galactic halo 

Annihilation of dark matter particles of 
mass m s and number density n{r) gener- 
ates a diffuse cosmic ray flux at our loca- 
tion r [31] 
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Here v = 1/2 if the annihilating parti- 
cles are Majorana to avoid overcounting 
of collisions [32], and v — 1/4 otherwise 



(assuming in the non- Majorana case equal 
amounts of dark matter and anti-matter 
in the halo). The case of interest to us is 
v = 1/2. The fragmentation function 



F(E, E in ) = 



<W(E, E m ) 
dE 



1 da 
adE 



gives the number of particles per energy 
interval and per annihilation event for 
an event with initial energy E in . Later 
on we will mostly focus on the photonic 
part dAf^/dE = a~ 1 da^ /dE, as charged 
and hadronic components of dark mat- 
ter annihilation products are masked by 
a relatively larger cosmic ray background 
than photons or neutrinos. However, in 
the present section we will also discuss the 
fraction of the total flux of cosmic rays 
from dark matter annihilation products 
compared to the cosmic rays flux jcR from 
all sources. 

The factor ov in equation (3) is the prod- 
uct of annihilation cross section and rela- 
tive speed in the non-relativistic limit. For 
processes at high redshift or annihilation 
of very low mass particles, thermal averag- 
ing (avj 7 ) would have to be included. For 
high redshift sources, redshifting of da /dE 
would also have to be included and the dis- 
tance \r & — r\ has to be replaced by the 
luminosity distance. However, for annihila- 
tion of heavy particles in the galactic halo, 
equation (3) and its corresponding line-of- 
sight counterparts below are perfectly ade- 
quate. 

Equation (3) yields a cosmic ray flux av- 
eraged over all directions. If the detector is 
only sensitive to cosmic rays from a small 
solid angle Af2 <C 1 sr, or if corresponding 
cuts can be applied, the observed flux per 
unit of solid angle is [33] (see also [34]), 
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The vector x with length x and direction 
p is related to the vector r in (3) through 
x = f — r* , n(x) = n(r). 

The averaged flux (3) is recovered from 
equation (4) in the following way. The ob- 
served flux for aperture AQ — > is the in- 
tegral along the line of sight ($, (p) , 



jAn-o(#, <p)= I dx 





vn (x, tp) 
Aix 



<W(E, 2m s ) 

x —— -av. 

dE 



(5) 



Averaging jAn->o(^ ; V 9 ) over an directions 
yields the diffuse flux (3), 



j = 0'An->o(A^)) 
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We will use equation (3) to estimate the 
averaged cosmic ray flux from dark matter 
annihilation in the galactic halo and equa- 
tion (5) to estimate the diffuse cosmic ray 
flux along a line of sight orthogonal to the 
galactic plane. 

We refer to the differential fluxes (3,4) as 
inclusive fluxes. Stable final states include 
photons, neutrinos, electrons, positrons, 
protons and anti-protons. The differential 
photon fluxes are found by substituting 
dU{E,2m s )/dE -> dN 1 {E,2m s )/dE in 
equations (3-5). 



The density profile is assumed as an 
NFW profile [35], 



p(r) 



M 



r(r + r s y 



(7) 



with a mass parameter M = 4.85 x 
1O 1O M = 5.41 x 10 64 TeV/c 2 and a scale 
radius r s = 21.5 kpc. These parameters 
correspond to the fit by Klypin et al. to 
the galactic halo [36], see also [37]. 

Equation (3) can be evaluated analyti- 
cally for an NFW profile (we use v = 1/2 
in the following), 



r + r \ p 2 (r) 




\ r — r e\ 

dAf(E, 2m s ) av 

x 

dE 47r sr 

M 2 dAf(E,2m s ) av 



4m| 




where Lo denotes Euler's dilogarithmic 
function^] [38,39]. 

Substitution of data for the galactic halo 
yields 



\14 



j = 1.95 x 10 

dU(E, 2m s 



TeV V 
m s J 
av 



x ■ 



dE cm 5 sr 
1.95 x 10 14 x (—J 



V m s 



2 Note that L 2 (x) = f(l-x) in [39]. Although 
of no practical relevance, we point out that 
the singularity at r s = r cancels between the 
last two terms in equation (8). 
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dj\f(x, 2ms) av 



dx 



TeV cm 5 sr ' 



(9) 



with the scaled energy variable x = E /ms- 
Looking only along a line of sight or- 
thogonal to the galactic plane to minimize 
background effects will cost approximately 
a factor 3 in observed flux from the galactic 
halo, 



is /TeVV 
6.20 x fO 13 x [ 



x 



V m s 
dj\f(x, 2ms) w 



dx 



TeV cm 5 sr 



(10) 



Cirelli et al. also compared their calcu- 
lations with other density profiles [20], us- 
ing the assumption of same local dark mat- 
ter density. This yields for the isothermal 
cored profile 



Q 



/i 



r 2 + r z r 



5 kpc, 



a parameter = 6.92 x 10 62 GeV/kpc = 
2.24 x 10 41 GeV/cm. 

The absence of the central cusp reduces 
the flux averaged over all directions, 
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8m|r c (r 2 + r 2 
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TeV cm 5 sr ' 



(11) 



but the flux along a line of sight perpendic- 
ular to the galactic plane is virtually un- 
changed and slightly increased due to the 
weaker local gradient in the dark matter 
distribution, 



jf c) = 6.36xl0 13 x /TeVV 



x 



V m s J 
dj\f(x, 2m s ) crv 



dx 



TeV cm 5 sr 



(12) 



We will report numerical results for the flux 
j from equation (9). The other fluxes can 
then easily be derived from equations (10- 
12). 

The fragmentation function into all final 
states must satisfy the energy sum rule 



/ 



dxx 



dM(x) 
dx 



An often used parametrization is 



dM(x) 



X) 



dx 



2x a (l 



B(a + 2,P + 1) 



(13) 



We will use fiducial values a = — 1.5, (3 = 2 
for numerical estimates. The two slope be- 
havior of \og(dM {x) / dx) versus x (see e.g. 
Fig. 17.1 in the review [40]) also suggests a 
phenomenological fit 



dAf(x) 
dx 



Aexp(-cxx) + Bexp(-f3x). (14) 



The values found by the HRS Collabo- 
ration in e + e~ annihilation at ^/s = 29 
GeV correspond to normalized values 
A = 360.5, a = 28.4, B = 97.2 and 
(3 = 7.91 [41]. A disadvantage of these 2- 
temperature distributions is that they are 
very small, but do not vanish at a; = 1. 
However, they work very well between 
0.1 < x < 0.9. 

The shape of fragmentation functions is 
only weakly energy dependent between 12 
GeV and 202 GeV [40], and dominated by 
the fragmentation properties of intermedi- 
ate partons. Therefore we also use (14) for 
numerical work besides (13). 

We will use the differential photon frag- 
mentation function proposed in [42], 



dJV^x) _ 0.42exp(-8. 



x 



dx 



x L5 + 0.00014' 



(15) 



for the differential photon spectrum. This 
corresponds to about 11% photon energy 
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count through a correction factor 



3 Model independent bounds on the 
cosmic ray flux from dark matter 
annihilation 

Although we are primarily interested in 
cosmic ray fluxes from the Higgs portal 
coupling (1), we would also like to point 
out that model independent estimates for 
cosmic rays from dark matter annihilation 
arise from unitarity limits and from limits 
on neutrino fluxes. Another model inde- 
pendent limit arises from halo stability 
[43]. However, this limit has been super- 
seded by the neutrino limit for dark matter 
masses heavier than 0.1 GeV [44,45]. 

The unitarity limit on s wave annihila- 
tion cross sections a < 4n/k 2 [46,47] im- 
plies 



av<AA0 x KT 19 



cm 



x 



/TeV\ 2 lOOkm/s 
V m s J v 



(16) 



Substitution in equation (9) yields a limit 
on the diffuse cosmic ray flux from galactic 
dark matter annihilation 



3s 



< 



8.57 x 1(T 5 dAf(x) 
Te V cm 2 s sr dx 
/TeV\ 5 lOOkm/s 

V m s J v 



:i7) 



Beacom et al. and Yiiksel et al. recently 
found that limits on the diffuse cosmic neu- 
trino signal and the halo signal can be used 
to impose stronger limits on the dark mat- 
ter annihilation cross section for dark mat- 
ter masses below 10 TeV [44,45]. Between 
10 GeV and 10 TeV, the cosmic and the 
isotropic halo neutrino flux limits approxi- 
mately reduce the limit on the annihilation 
cross section ~ m 2 . We take this into ac- 



(™) 2 > 10GeV <™s<10TeV, 



m s > 10 TeV. 



Comparison with the cosmic ray flux be- 
low 1 PeV [48] 



3CR ■ 



2.582 x 10" 5 /TeV\ 
TeV cm 2 s sr \ E ) 



2.68 



2.68 



_ 1.236 x 10~ 5 / 100 GeV \ 
GeV cm 2 s sr V E J 

shows that a galactic dark matter annihi- 
lation signal could reach several per cent of 
the total differential cosmic ray flux if the 
annihilation cross section could get close to 
the upper limits either through Sommer- 
feld enhancement or through resonance ef- 
fects. We find 



3s 

3CR 



< A, x 3.32a; 



2.68 



d/V(x) 
dx 



/ TeV \ 232 lOOkm/s 
V m s J v 



e.g. for m s = 100 GeV, 



3s 
3cr 



< 6.94 x 10 - V- 68 



d/V(x) lOOkm/s 



dx 



x 



The maximum of x 2S8 d/V(x)/dx for the 
fragmentation function (13) is 0.23 for x = 
0.37. For the fragmentation function (14) 
the maximum of x 2m dj\f{x) / dx is 0.37 for 
x = 0.34. 

However, cosmic rays are strongly dom- 
inated by charged particles, while dark 
matter annihilation products are expected 
to contain a relatively higher neutral com- 
ponent of photons and neutrinos. 7 ray and 
neutrino observatories are therefore the 
primary search tools for dark matter anni- 
hilation products. We will use the diffuse 
photon background published by EGRET 
between 30 MeV and 120 GeV [49], 
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2.10 



. _ 6.89 x 1(T 10 /TeVV 
HE ~ TeVcm 2 ssr V~zT/ 

8.68 x 10- 11 /100GeV\ 210 



GeV cm 2 s sr 



E 



for benchmarking. This is the diffuse pho- 
ton flux observed by EGRET after subtrac- 
tion of conventional galactic sources. It has 
been pointed out that improved models for 
standard sources could reduce the EGRET 
signal, see [50] and references there. The 
sensitivity calibration beyond 1 GeV has 
also been called into question [51]. The en- 
ergy range between 1 GeV and 100 GeV 
is particularly relevant for continuous pho- 
ton signals from dark matter in the mass 
range of interest here, 80 GeV < ms < 
1 TeV, and we will see that signal levels 
may be small. Reliable further subtractions 
of standard sources from the diffuse 7 ray 
background signal, and a lower total sig- 
nal, would help to identify or constrain a 
possible dark matter signal. EGRET sets a 
useful benchmark until GLAST/LAT pub- 
lishes data on the diffuse 7 ray background. 

Normalizing the averaged photon flux 
following from equation (9) to the EGRET 
flux (18) yields 



^L = 2.83x loV 1 ^^ 
He ax 



V m s ) 



(TV 



cm 3 s 1 



4 Annihilation cross sections for the 
electroweak singlet 

The coupling (1) reduces in unitary 
gauge to 

Hsh = T±j d z SS 2 h + || d?xS 2 h 2 . (20) 

For the annihilation of the electroweak sin- 
glet through an intermediate Higgs, we also 
need the couplings 



H, 



hh 



,^mt / o h 





H fh = fd 3 xJ2^hf.f 



and 

Hw,zh 



2^W'W + + 



(21) 
(22) 



v h 



x K^ + y 1 ■ 



(23) 



These couplings yield the following 
annihilation cross sections in the non- 
relativistic limit, 



vass^hh = V 



x 



W - rn h 2 
16"7rm,s 3 

2m s 2 + m h 2 



4m s 2 - m h 2 + im h r h 



2r]v h 2 



2m s 2 — m h 2 



The maximum of x 2 - 1 dAf 1 (x)/dx is 0.05 
near x ~ 0.08. Therefore any annihilation 
cross section of a dark matter particle with 
a mass below 1.5 TeV is constrained to 



va 



ss^ff - V 



2 N c m f 2 

4:71171s 3 



X 



ms 2 — rrif 2 



(Am s 2 - m h 2 ) 2 + m h 2 T h 2 ' 



(TV < 10" 



-2:-! 



cm 



x 



UOOGeV/ 



0.9 



(19) 



va ss ^ww = V 



, \/m s 2 — mw 2 
4nm s 3 
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3m w 4 — Am w 2 m s 2 + Ams A 
(Am s 2 - m h 2 ) 2 + m h 2 T h 2 



vcrss^zz = V 



,\fm s 2 — m z 2 

8nm s 3 
3mz 4 — Amz 2 ms 2 + Ams A 
(Am s 2 - m h 2 ) 2 + m h 2 T h 2 ' 



with N c = 3 for quarks and N c = 1 
for leptons. The & S s^fj cross section is 
summed over final spin states. The cross 
section ass^hh contains the scattering am- 
plitude from the S 2 h 2 contact vertex in 
equation (20), the scattering amplitude 
from the s-channel contribution with an 
intermediate Higgs boson from the S 2 h 
vertex in equation (20) and the h 3 vertex 
in equation (21), and the t-channel and 
li-channel amplitudes with an intermedi- 
ate singlet from the S 2 h vertex. Those 
amplitudes are with the normalization 
S f i = 5 fi - iM fi S A (p 1 +p 2 -h- h), 



M {1) - V 



16tt 2 

X 



^EsfyEsfyEhip^Ehipz) ' 
3r]m h 2 



lQn 2 



x 



x 



M SS-*hh ~~ 8n 2 



E s {k 1 )E s {k 2 )E h {p 1 )E h {p 2 ) 
1 

(ki + k 2 ) 2 + m h 2 - ie' 

7] 2 V h 2 



X 



X 



+ 



\J Es(ki)E s (k 2 )E h (pi)E h (p 2 ) 

( I 

\(^i ~ Pi) 2 + ms 2 - ie 

! ). 

(h -p 2 ) 2 + m s 2 - ie/ 



Note that due to the constraint mj 2 > 
mh 2 for SS — > hh, neither the s-channel 
resonance for ms = mh/2 nor the t-channel 
or -u-channel resonance for ms = rrih/^/2 
can be realized for oss^hh- Furthermore, 
Ts = in our models. 

The amplitude for annihilation of two 

— * — * 

singlets with momenta k\ and k 2 into mas- 
sive vector bosons with momenta p±, p 2 
and polarizations ^(pi), e^(p 2 ) is pro- 
portional to mw 2 , 



M 



ss^ww 



X 



X 



8tt 2 
1 



Esik^Esik^Ewip^Ewip,) 
e M(^). 6 (/3)(^ 2 ) 

(ki + k 2 ) 2 + m h 2 - ie 



(24) 



However, the longitudinal polarization vec- 
tor for massive vector bosons comes with a 
factor mw" 1 , e.g. the polarization vectors 
for p in ^-direction are 

e<»(p) = (0,1,0,0), e<> 2 \p) = (0,0,1,0), 
e (3) (p1 1 



m w 



1^,0,0,^ + ^2 



The tensor product of massive polarization 
vectors 



J2^(p)^^ a)u (p) = rr + 

implies 



a,/3 

\ m w I V m w z 



= 2 + 



(Pi ■ P2? 

m w 4 



Together with the amplitude (24) this 
yields the cross section ass^ww- An equiv- 
alent way to understand why the cross 
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sections a S s^ww and &ss->zz do not 
vanish in the limit of vanishing coupling 
Vh 2 ~ my/ 2 ~ mz 2 — > are the residual 
Goldstone boson modes /i 1 * 1 and x° from 
the Higgs field, which would contribute to 
singlet annihilation in the SU(2) x U(l) 
symmetric limit. 

Standard model decay widths of the 
Higgs particle are small unless the Higgs 
boson is heavy enough to decay into 
weak gauge bosons [52], e.g. Yhijrih = 
115 GeV) = 4.56 MeV, T h (m h = 160 GeV) = 
6.37 MeV, T h (m h = 165 GeV) = 247 MeV, 
r h (m h = 203 GeV) = 1.54 GeV. 

The energy dependent cross sections 
av(K) in order 0(rf) in all cases corre- 
spond to the substitution 

m 2 s -> (K + m s f 

in the equations for va, where K is the 
kinetic energy of an electroweak singlet. 
The non-relativistic limits for va are there- 
fore not only excellent approximations for 
T <C ms, but also provide upper limits on 
the thermal averages at high temperature 



(av) 



Ills i 

JO alX cxp[(K+m s )/T\-l 



The temperature dependence is very 
weak. The 0(rf) cross section for m s = 
200 GeV is only reduced from av = 
1.04?7 2 x 10" 23 cm 3 /s at T = to (av) = 
0.39r] 2 x 10~ 23 cm 3 /s at T = 100 GeV. 

The tree-level cross sections reported 
here yield cross sections of order av ~ 
rf x 10~ 23 cm 3 /s in the mass range 
m s ~ m h , see also figures 1-5 below. In ad- 
dition there are also loop suppressed pho- 
ton lines at E 1 = ms from a 77 final state, 
and at E 1 = m s — (mz 2 /<ims) from a 7Z 
state. These contributions are of order of a 
per cent compared to the continuous spec- 
trum, because the loop amplitudes contain 



two vertices of electroweak strength or the 
direct W + W~ — > 77 contact vertex. The 
contribution to v ass->j-y through a W + W~ 
loop with contact vertex e.g. is approx- 
imately of order (neglecting logarithmic 
mass dependencies) 



va S s- 



32ira i] mw 



>77 m s 2 [(4m 5 2 - m h 2 ) 2 + m h 2 T h 2 ] 



10- s va SS - 



The continuous spectrum yields a con- 
servatively estimated photon energy yield 
of order 10% (see equation (15)). Together 
with the ratio of cross sections, this im- 
plies that in terms of photon energy yields 
from singlet annihilation, the line spectrum 
should be of order of a few per cent of the 
continuous spectrum. Therefore we focus 
on the continuous spectrum in the present 
investigation. However, further study of the 
loop induced line spectrum is also of inter- 
est. 



4-1 Annihilation signal from a heavy elec- 
troweak singlet in tree level approxi- 
mation 

If we assume a Higgs mass limit < 
203 GeV from electroweak analysis [53] 
and absence of a very heavy fourth gen- 
eration below ms, we get in leading order 
the following annihilation cross section for 
a heavy (ms > 1 TeV) electroweak singlet 
coupling through the Higgs portal, 



va S s^ 



7 V 2 



64:7im s 
2 



25 /TeV\ 2 cnr 



= 4.06?T x 10"^ 

V ms J 



■ (25) 



This translates into a cosmic ray flux 
from heavy electroweak singlet annihila- 
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tion through the Higgs portal at tree level, 

. _ 7.92r/ 2 x 10~ n /TeV\ 5 dM{x) 
^ S TeVcm 2 ssr Vm^j dx 

This is very small compared to the galac- 
tic background flux from supernovae and 
their remnants, and again photon signals 
are expected to be more sensitive. Normal- 
izing to the EGRET flux formally yields 

A. = o.ii5,v-'^M (™)". 

h,E dx Us / 

Due to the energy limit of EGRET this 
should only be used for x < 0.12 TeV/ms- 
For ms = 1 TeV and E ~ 76 GeV we find 
h/h,E ^ 5.6/7 2 x 10~ 3 . 

For heavier singlets, March-Russell et 
al. have pointed out that Sommerfeld en- 
hancement due to scalar exchange between 
annihilating supersymmetric singlets can 
boost annihilation cross sections by factors 
10 3 to 10 5 for singlet masses between 1 TeV 
and 30 TeV [30] . Sommerfeld enhanced an- 
nihilation signals from heavy electroweak 
multiplets have recently been discussed 
by Cirelli et al. [20]. In the present paper, 
we will instead focus on electroweak sin- 
glet annihilation in the intermediate mass 
range 80 GeV < m s < 1 TeV, when ass 
is enhanced due to proximity to the W, Z 
and Higgs peaks in the cross section. 

4-2 Enhancement of electroweak singlet 
annihilation for ms ~ rrih due to the 
W , Z and Higgs peaks 

Electroweak singlet annihilation through 
an intermediate Higgs can be strongly en- 
hanced when the channels SS — > WW and 
SS — > ZZ open up. This is especially rele- 
vant when mj ~ rrih, which is the primary 
mass range of interest in our present inves- 
tigation. 



We will primarily use Higgs mass val- 
ues rrih = 115 GeV from the direct search 
limit [53], and rrih = 160 GeV in the range 
of highest sensitivity for search for a light 
Standard Model Higgs boson at the Teva- 
tron [54]. The value rrih — 160 GeV is also 
in the preferred mass range for minimal 
dark matter models identified in reference 
[14]. 

Compared to ass^ww the contributions 
of light quarks and leptons to a S s are of 
order 10~ 4 and the contributions of heavy 
quarks and leptons are of order 10 -2 in the 
mass range of interest. We will include c, b 
and t quarks with masses m c = 1.25 GeV, 
m h = 4.20 GeV and m t = 172.5 GeV in 
the calculation of v ass- The r lepton will 
be included with m T = 1.78 GeV. 
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80 90 100 110 120 130 140 

m s JGeV 



Fig. 1. The cross section n~ 2 voss * n units of 
10 -23 cm 3 /s for very weak coupling rf < 0.01. 
The red (initially and finally highest) curve is 
for nih = 129 GeV, the blue (initially and fi- 
nally middle) curve is for rrih = 115 GeV, and 
the green (initially and finally lowest) curve is 
for m h = 89 GeV. 

The effect of the t-channel plus -u-channel 

(3) 

amplitude M-ss-^hh ^ s sman f° r ver Y weak 
coupling rf < 0.01. The Oijf) cross sec- 
tions for very weak coupling are displayed 
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for several Higgs mass values in figures 1 
and 2. 

Figure 1 shows the scaled cross sec- 
tion r]~ 2 vass for = 115 GeV and 
70 GeV < m s < 140 GeV. The cross sec- 
tions for rrih = 89 GeV and m^ = 129 GeV 
are included for comparison. The edge at 
80 GeV arises from the opening of the W 
channel SS -> WW, the edge at 91 GeV 
arises from the Z channel, and the edge at 
ms = rrih arises from the Higgs channel 
SS -> hh. 
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120 130 140 150 160 170 180 190 200 

m^GeV 



Fig. 2. The cross section r]~ 2 vass in units of 
10 -23 cm 3 /s for very weak coupling rj 2 < 0.01. 
The blue (initially and finally lower) curve is 
the cross section for = 160 GeV. The red 
(initially and finally upper) curve is the cross 
section for m^ = 165 GeV. 

The scaled cross sections r]~ 2 v<Jss for 
rrih = 160 GeV and rrih = 165 GeV are 
displayed in figure 2. 

The cross sections for larger values of ms 
approach the asymptotic limit (25), see fig- 
ure 3. 

For electroweak strength coupling rf = 
0.1, the t-channel plus it-channel amplitude 
M-ss-*hh suppresses the Higgs threshold in 
the cross sections, see figures 4 and 5. 

The continuous photon spectrum for 
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Fig. 3. The cross section n~ 2 vass in 
units of 10 _23 cm 3 /s in the mass range 
160 GeV < ms < 640 GeV for very weak cou- 
pling rj 2 < 0.01. The red (initially upper) curve 
is the asymptotic cross section (25). The blue 
(initially lower) curve is for = 115 GeV. 
The cross section for rrih = 160 GeV also dif- 
fers by less than 10% from the aymptotic limit 
for m s > 400 GeV. 

£ 7 < 120 GeV in units of the EGRET 
flux is displayed in figure 6 for mj = 120 
GeV, m h = 115 GeV, rj 2 = 0.1 and 
vo- ss = 2.03 x 10" 24 cm 3 /s. 

The corresponding flux ratio for ms = 
120 GeV, m h = 160 GeV, rj 2 = 0.1 and 
v&ss = 3.87 x 10~ 24 cm 3 /s is diplayed in 
figure 7. 

The fluxes in figures 6 and 7 assume that 
all cold dark matter is in a singlet state, 
but do not take into account correlations 
between ms and mh from thermal creation 
of singlets. This topic will be addressed in 
the next section. 

5 Thermal creation and abundance 
estimates 

A thorough and beautiful application of 
Lee- Weinberg theory [55] for abundance 
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70 80 90 100 110 120 130 140 
m s /GeV 

Fig. 4. The cross section vo~ss ^ n units of 
10 -24 cm 3 /s for coupling rf = 0.1. The red 
(initially and finally highest) curve is for 
nth = 129 GeV, the blue (initially and finally 
middle) curve is for = 115 GeV, and the 
green (initially and finally lowest) curve is for 
m h = 89 GeV. 

estimates of electroweak singlets has been 
given in [12]. We will revisit the subject for 
the particular range of masses and cross 
sections of interest to us. The applica- 
tion of Lee- Weinberg theory in reference 
[12] used Vtdm = 1 while the application 
for light electroweak singlets in [13] used 
fldm = 0.6. Here we use VL dm = 0.2 [53]. 

In the absence of more detailed model 
assumptions about singlet generating in- 
teractions behind the Higgs portal or co- 
annihilations, the effect of thermal cre- 
ation of electroweak singlets below the 
electroweak phase transition is taken into 
account through a thermal production 
term in the rate equation 

|(na 3 ) = Nthermal ~ (<™>nV , (26) 

where a(t) is the scale factor in the 
Robert son- Walker metric. The thermal 
production term is determined from the 



120 130 140 150 160 170 180 190 200 

m s /GeV 

Fig. 5. The cross section vass i n units of 
10 _ cm 3 /s for coupling rf = 0.1. The blue 
(initially and finally lower) curve is the cross 
section for nih = 160 GeV. The red (initially 
and finally upper) curve is the cross section for 
m h = 165 GeV. 




oH ■ 1 ■ 1 ■ 1 i — — i 

0.2 0.4 0.6 0.8 1 

E ,/ m s 

Fig. 6. The photon flux j 7 in the continuous 
spectrum for mh = 115GeV ; ms = 120 GeV, 
if = 0.1 and E~ < 120 GeV in units of the 
EGRET flux. 

equilibrium requirement d(na 3 )/dt = 
[55,56], 

Nthermal = {(Tv)^^ , (27) 



12 



0.3 ■ 



0.2 



0.1 



0.2 0.4 0.6 

E I m,. 

7 J 



Fig. 7. The photon flux j 7 in the continuous 
spectrum for nth = 160 GeV, ms = 120 GeV, 
rj 1 = 0.1 and E y < 120 GeV in units of the 
EGRET flux. 

where Uq is the thermal equilibrium den- 
sity for temperature T. The resulting rate 
equations during radiation domination is 

dn 3n , . / , 9 \ , . 

* + 2F = - (OT 'X» 2 -«»)' < 28 > 

Radiation domination also yields 
b 



t 



Ji2 ' 



(29) 



and this is used to rewrite equation (28) in 
the form 



d n , . . n 2 — nk 



(30) 



Equation (29) follows from the relations 
for the energy density in radiation for 
Unfiation <t < t eg (t eq ~ 2.4 x 10 12 s is the 
time of radiation-matter equality), 



£ 7 = 9*{ T ) 



7T 



'{ksT) 4 3mp lanck c 



30(hc) 



4ht 2 



(31) 



Here we use the reduced Planck mass 
mpi anck = (hc/87rG N ) 1/2 . We have g*(T) = 
91.5 for 777.5 < T < mw The parameter b is 



ShmpianckC 2 



nk 2 B 



= 2.53 x KT'fci^sGeV 
= 3.41 x 10 19 sK 2 . 

The analytic approximation proposed by 
Lee and Weinberg [55] uses the equilibrium 
density 



n (T) 



1 



27T 2 (^C) 



7 (K + mc 2 )JK(K + 2mc 2 ) 
x / dK- 1 

o 



exp[(K + mc 2 )/k B T] - 1 



until a freeze-out temperature Tf is reached 
with 



d n (T) 



dT T 3 



2b(av) 



n 2 (T f ) 



(32) 



This temperature usually turns out to 
satisfy Tf < 0.05ms, such that the non- 
relativistic limit of no(T), 



n (T) 




m s k B T 



2tt 



exp 



rrisc 
~T^T 



can be used in the evaluation of the Lee- 
Weinberg condition (32). The solution is 
then extended for T < Tf using domina- 
tion of the annihilation term in the rate 
equation (26), 



dn 3n 
~dl + 1 2l 



-(av)n 2 , t f <t<t eq (33) 



dn 3nda , . ~ 
— + — — = -(crv)n 2 , t > t, 
dt a dt 



(:'(/■ 



(34) 



with the initial condition n(tf) = n (Tf). 
Here a = a(t) is the scale factor in the 
Robertson- Walker line element. One might 
use the following set of equations for t > 

teq j 

dn 2n 2 
— + — = -{av)n , t eq <t < t A , 



dt 



t 
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dn 2n _ / t \ . . 2 

— H coth — = —{av)n , t > t A , 

at r A Vta/ 

with the time constant (for A = 0.76g c = 
4.27keV/cm 3 ) 



ta — 2mpi anc k, l , 
V 3nA 

and £ A following from 



3.23 x 10 17 s, 



1 + ^ 



a(f ) /sinh(t /T- A ) 



2/3 



a(t A ) \smh(t A /T A ) 



/Qa 



1/3 



1.47, 



see e.g. the appendix to reference [57] . How- 
ever, we will see from the solution of equa- 
tion (33) that n(t eq )(av)t eq < 1(T 9 < 1, 
see equation (39) below. Therefore both 
n(t){av)t <C 1 andn(t)(<Tf )rAtanh(t/rA) <C 
1 for t > t eq , i.e. the expansion term dom- 
inates strongly over the annihilation term 
for t > t eq . This yields standard cold dark 
matter evolution for late times, 



n(t) = n(t 



' a(t 



eq ) 



eq) 



a(t) 



t > t 



eq- 



We define £ = mgc 2 /ksTf. The Lee- 
Weinberg condition (32) takes the follow- 
ing form 

2bk 2 B m s c 2 
exp(0= ; 7^=r \ av ) 



(v^vrftc) 3 
4.18 x 10 11 



£-1.5 

(av) 



x ■ 



rrtsc 



10 -24 cm 3/ s 



100 GeV£- 1.5 



(35) 



We are interested in a perturbative 
Higgs portal rf < 0.1 and weak scale 
singlet masses. Our previous results on 
cross sections then imply that the fac- 
tor (((Tu)/10- 24 cm 3 s- 1 ) x (m s /100GeV) 
should be in the range between 0.1 and 10. 



This will yield values for £ between 20 and 
30. Thermal theories of particle creation 
using equation (27) generically predict 
that particles will remain thermal until the 
temperature has dropped to a value well 
below their mass threshold. 

Integration of equation (33) yields 



n(t 



eq) 




l + 2n(t f ){av)t f \t eq J 



(36) 



where we used that freeze-out tempera- 
tures following from (35) will at least be a 
few GeV for the parameter range of inter- 
est here, and therefore {t eq jtf) x l 2 > 10 9 . 

The relation between temperature and 
time and the definition of £ imply 



bk B 2 , 2 



Tf 2 



ms 2 c A 



= 2.53 x 10^ n £ 2 s x 



100GeV\ 2 



m s c* 



(37) 



This yields 



(<ru)f / = 2.53 x 10~ 35 £ 2 cm 

/100GeV\ 2 

x —) x 

V m s c z J 



(av) 



10 -24 cm 3/ s - 
On the other hand, the density is with (35) 

/ \ 3 

m s c 



n 



(*/) = 



V2^h / 

1 4 
Ills C 



exp(-£) 



1.5 



2bk B 2 (av) £ 2 
1.98 x 10 34 



£-1.5 



ra$c 



e 

2 



100 GeV 



x 



-cm 3 

1Q -24 



cm 3 /s 



(av) 
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and therefore 
2n(tf){av)t f = S -1.5. 
This equation also implies 

n(t f )(av)tf 



(38) 



n(t eq )(av)t, 



eq - 



l + 2n(tf)(av)t f V t 




£-1-5 
2£-l 



'^<io- 9 , 



(39) 



and therefore the annihilation term is much 
smaller than the expansion term for t >t eq , 
d(na 3 )/dt = 0, 



n(to) = n{t eq )z eq 3 , 



(40) 



where 



a(t ) 

a{teq) 



a(to) 
a{t eq ) 



We use z eq = 3000 for numerical work. 

Equations (36), (37), (38), and (40) yield 
the current energy density in one singlet 
species, 



Qs ] =n{t eq )m s c 2 



k B Vb 



2i-r2{av)t eq ^z eq z 
2^-3^ 2.51eV/cm 3 



2£-r (av}/10- M cm 3 /s' 



(41) 



However, we can have an O(N) symmetric 
iV-plet of electroweak singlets of mass m S) 
such that the current energy density in sin- 
glets is 



Qs 



2£ 



3, iV x 2.51eV/cm 3 
< x — /,n-^ --.,/.■ • (42) 



2£-r (av)/10- 24 cm 3 /s' 



The condition g s = g dm = 1.106 keV/cm 3 
then determines £ in terms of (av)/N, 
and substitution in equation (35) then 
relates m s , rrih, and N. We report the re- 
sulting singlet masses in the mass range 





m h = 


= 115 GeV 


m h = 


= 160 GeV 


N = 


1 


m s = 


= 879 GeV 


m s = 


= 883 GeV 






Tf = 


35.9 GeV 


Tf = 


36.1 GeV 






X(j — 


0.0531 


x a — 


0.0531 


N = 


10 


m s = 


= 273 GeV 


m s = 


= 285 GeV 






Tf = 


10.7 GeV 


Tf = 


11.1 GeV 








0.557 




0.558 



Table 1: Singlet masses in the mass range 
mw < ms < 1 TeV which satisfy the Lee- 
Weinberg condition (35) for rf = 0.1, = 
115 GeV or m h = 160 GeV, and N = 1 or 
N = 10. 





m h = 


- 115 GeV 


m h = 


= 160 GeV 


N = 1 


m s = 


= 293 GeV 


m s = 


= 304 GeV 




Tf = 


12.6 GeV 


Tf = 


13.0 GeV 






0.0505 


X(j 


0.0506 


N = 10 


n/a 




m s = 


= 111 GeV 








Tf = 


4.47 GeV 








x a — 


0.540 



Table 2: Singlet masses in the mass range 
mw < ms < 1 TeV which satisfy the Lee- 
Weinberg condition (35) for rj 2 = 0.01, = 
115 GeV or m h = 160 GeV, and N = 1 or 
N = 10. 

m w < m s < 1 TeV for rf = 0.1 and for 
iV = loriV = 10 in table 1. We also report 
the corresponding freeze out temperatures. 
The annihilation cross sections are given 
in the form 

vo- ss 



10-24 cm 3/ s - 

Table 2 shows solutions in the mass range 
m w < m s < 1 TeV for very weak coupling 
rj 2 = 0.01. 

The cross sections in table 2 for given 
mh and iV are similar to the cross sections 
in table 1, in spite of the weaker coupling. 
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This is due to the fact that the correspond- 
ing singlet masses in table 2 are smaller and 
much closer to the peaks in the cross sec- 
tion. 

For the flux calculations for N > 1, we 
have to rescale the flux (9) by a factor 1/N, 
because the density factor n 2 for each anni- 
hilating species is now suppressed oc iV~ 2 , 
but there are N annihilating species of sin- 
glets of mass ms- Therefore we find with 
equations (9) and (15) 



1.95 x 10 20 /GeV\ 

x 

N Us 

dM 7 (E y ) va ss 



■ 



cur sr 



8.19 x 10 19 /GeV\ 



x 



') 



X 



N V ms 

v /m^exp(-8E 7 /m s ) va ss 
£L L5 + 0.00014m s L5 cm 5 sr' 



(43) 



The parameter ( vass /./V) / ( 1 0~ 24 cm 3 /s) - 
x a /N varies only in the range 5.05 x 10~ 2 < 
x a /N < 5.58 x 10~ 2 for the solutions in 
tables 1 and 2. However, j 7 scales with 
m^" 1,5 , and therefore the low-mass solu- 
tions from table 2 yield higher flux than 
the solutions from table 1. We will give 
results for j 7 in the case rf = 0.1, iV = 10, 
m h = 160 GeV and m s = 285 GeV, 
and also in the case rf = 0.01, N = 10, 
m h = 160 GeV and m s = 111 GeV. 

The expected contribution to the photon 
flux below 120 GeV in units of the EGRET 
flux is shown in figures 8 and 9 for the two 
cases. The corresponding number of pho- 
tons per GeV • cm 2 • s • sr is given in figures 
10 and 11, respectivley. 

The integrated photon flux 



$ 7 (£ 7 ) = J dE h (E) 



0.0020 



0.0015 



0.0010 



0.0005 



0.1 
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E Irric 

T * 
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0.4 



Fig. 8. The photon flux j 7 in the contin- 
uous spectrum for rf = 0.1, N = 10, 
m h = 160 GeV, m s = 285 GeV and 
£ 7 < 120 GeV in units of the EGRET flux. 



0.005 



0.004 



0.003 



0.002 



0.001 




in the two cases is shown in figures 12 and 



Fig. 9. The photon flux j 7 in the contin- 
uous spectrum for n 2 = 0.01, N = 10, 
m h = 160 GeV, m s = 111 GeV and 
E~ t < 111 GeV in units of the EGRET flux. 

13. 

For an instrument with an effective area 
of 8000 cm 2 and a field of view of 2.4 sr, 
comparable to the Large Area Telescope 
on GLAST, the photon flux j 7 in figure 
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Fig. 10. The photon flux j 7 in units of 
(GeVcn^ssr)^ 1 for r/ 2 = 0.1, N = 10, 
m h = 160 GeV, m s = 285 GeV and £ 7 be- 
tween 15GeV and 30 GeV. The maximum in 
figure 8 corresponds to EL ~ 22 GeV. 



2.2 x 10" 10 



2. x 10" 



1.8 x 10" 



1.6 x 10" 



1.4 x 10" 



1.2 x 10 



-10 



10 



1.x 10" 
8. x 10" 11 
6.x 10" 11 
4.x 10" 11 



10 
£/GeV 



12 



14 



Fig. 11. The photon flux j~ in units of 
(GeVcm 2 ssr)- 1 for rf = 0.01, N = 10, 
nih = 160 GeV, ms = 111 GeV and EL be- 
tween 5 GeV and 15 GeV. The maximum in 
figure 9 corresponds to E 7 — 9 GeV. 

11 corresponds to an annual rate of 470 
photons with energies between 5 GeV < 
EL < 15 GeV from dark matter annihila- 



Fig. 12. The integrated photon flux $ 7 above 



energy E 7 in units of (cm s sr) 
N = 10, m h = 160 GeV, m s 



i j or rf = q.1 ; 
= 285 GeV and 



E 7 between 100 MeV and 10 GeV. 
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Fig. 13. The integrated photon flux <3? 7 
above energy E~ in units of (cm 2 ssr)~ 1 for 



= 0.01, N ■- 
ms = 111 GeV anc 
10 GeV. 



10, m h 



160 GeV, 



E~ between 100 MeV and 



tion on a diffuse astrophysical background 
of 90,000 photons per year. The flux in fig- 
ure 10 corresponds to an annual rate of 48 
photons with energies between 15 GeV < 
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E 1 < 30 GeV from dark matter annihila- 
tion on a diffuse astrophysical background 
of 22500 photons per year. There are two 
characteristic features which would help to 
identify j 1 as an excess due to dark matter 
annihilation. The excess would extend over 
an energy range of order 10 GeV, and it 
would be correlated with the galactic halo. 

6 Conclusions 

We have considered perturbatively cou- 
pled electroweak singlet dark matter in 
the intermediate mass range my/ < ms < 
1 TeV. The leading order annihilation 
cross section of the singlets is enhanced 
and varies strongly due to proximity to 
the W, Z and Higgs peaks. The product 
vass is of order rf x 10 _23 cm 3 /s for mj 
close to m h , i.e. it can substantially exceed 
standard estimates of dark matter annihi- 
lation cross sections even for perturbative 
singlet-Higgs coupling. 

For singlet masses above the SS — > WW 
threshold and electroweak strength cou- 
pling rf ~ 0.1, the Lee- Weinberg condition 
and the requirement Qs = ^dm push mj 
to high mass values around 900 GeV if 
there is only one singlet state. However, if 
there is an iV-plet of electroweak singlets 
or if the coupling is weaker, rf < 0.01, 
lower singlet mass values can be achieved, 
and the annihilation signal from the con- 
tinuous 7 ray spectrum can reach a level of 
several per mil of the EGRET diffuse 7 ray 
signal for photon energies E 1 ~ 0.08m,g. 
This excess contribution over the expected 
cosmological background would appear 
typically over an energy range of order 10 
GeV, and its dark matter signature would 
be its correlation with the galactic halo. 

The flux reported by EGRET was the 
diffuse 7 ray flux after subtraction of then 
known or expected galactic components. 



Reduction of the diffuse "excess" 7 ray 
flux due to subtraction of a larger com- 
ponent from interstellar gas and standard 
extragalactic sources increases the relative 
importance of an annihilation signal for 
any possible excess signal and improves 
detectability. The diffuse 7 ray flux will 
be measured with higher sensitivity and 
precision in the near future by the Large 
Area Telescope aboard the GLAST satel- 
lite. This will also cover a larger energy 
range up to 300 GeV. 

The minimal dark matter models con- 
sidered here include four basic parameters, 
the singlet mass ms, the number of singlet 
states N, the singlet-Higgs coupling 77, and 
the Higgs mass m h . The assumption that 
electroweak singlets of mass m s provide the 
dark matter in the universe relates these 
parameters through Lee- Weinberg theory. 
Measuring the Higgs mass at the Tevatron 
or the LHC will reduce the number of free 
parameters in this class of minimal dark 
matter models to two, or maybe even to 
only one free parameter if a missing energy 
signal can be used to constrain a combina- 
tion of 77, N and m s . The anticipated small- 
ness of singlet annihilation signals in cos- 
mic 7 rays indicates that observation of the 
Higgs particle at the Tevatron or the LHC 
may be needed for a successful search for a 
singlet annihilation signal in cosmic rays. 
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